PROCEEDINGS 

OF SCIENCE 



Hi 

(N 



43 



in 



BONN-TH-2010-14, DESY 11-006, FR-PHENO-20 10-044 

Introducing HiggsBounds 2.0.0 



Philip Bechtle 

DESY, Notkestrasse 85, 22607 Hamburg, Germany 

Oliver Brein 

^v^j ' Physikalisches Institut, Albert-Ludwigs-Universitat Freiburg, 

Hermann-Herder-Str. 3, D-79104 Freiburg im Breisgau, Germany 



Sven Heinemeyer 

Institute de Fisica de Cantabria ( CSIC-UC), Santander, Spain 

Georg Weiglein 

43 

DESY, Notkestrasse 85, 22607 Hamburg, Germany 

Karina Williams 

Bethe Center for Theoretical Physics, Physikalisches Institut der Universitdt Bonn, 



Nussallee 12, 53115 Bonn, Germany 



J> ' We introduce version 2.0.0 of the computer program HiggsBounds! The program tests neutral 

O ' 



and charged Higgs sectors of arbitrary models against the current exclusion bounds from LEP 
and the Tevatron. As input, it requires a selection of model predictions, such as Higgs masses, 
branching ratios, effective couplings and total decay widths. The program uses the expected and 
observed topological cross section limits from the Higgs searches to determine whether a given 
parameter scenario of a model is excluded at the 95% C.L. by those searches. Version 2.0.0 
includes 39/53 LEP/Tevatron Higgs search analyses, compared to 13/36 in the previous release 
(1.2.0). Among the newly included analyses are LEP searches for neutral Higgs bosons (H) 
' decaying invisibly or into hadrons, LEP searches via the production modes t + t~H and bbH, 

Tevatron searches via ttH, and LEP and Tevatron searches for charged Higgs bosons. Also, all 
Tevatron results presented at the ICHEP'10 are included in version 2.0.0. In this note, we ex- 
plain the basic ideas behind the implementation of HiggsBounds and provide a list of search 
topologies implemented in version 2.0.0. Furthermore, we apply HiggsBounds 2.0.0 to (a) 
determine the allowed Higgs mass range for a simple 4th generation model, (b) update/reproduce 
LEP/Tevatron Higgs exclusion plots for the MSSM m™ ax benchmark scenario, and (c) show ex- 
clusion results for the scalar sector of the Randall-Sundrum model. 

Third International Workshop on Prospects for Charged Higgs Discovery at Colliders - CHARGED2010, 
September 27-30, 2010 
Uppsala Sweden 
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^For code download, online version, updated manual see : www . ippp . dur . ac . uk/HiggsBounds/ 
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1. Motivation 



The search for Higgs bosons, neutral or charged, is a major part in the endeavour to unravel the 
nature of electroweak symmetry breaking. The discovery of a charged Higgs boson, in particular, 
would be an unambiguous sign of an extended Higgs sector. The LEP experiments have searched 
for the Standard Model (SM) Higgs boson [jl|] and many others, and so do the Tevatron and LHC 
experiments. So far, no signals of Higgs bosons have been found, and LEP & Tevatron experimen- 
talists turned this non-observation into cross section constraints. The constraints are provided in 
the form of limits on cross sections of individual signal topologies (such as e + e~ — > h{Z — > bbZ or 
pp — > hiZ — > bbl + l~) or in the form of combined limits for a specific model, such as the SM. In the 
latter case, the individual topological cross sections have been combined assuming the proportions 
of the individual contributions to be as predicted by the model. 

The claim of model exclusion using non-observation constraints is a statement about the sta- 
tistical rejection of the "Higgs hypothesis" (i.e. that Higgs signal + background describe the data). 
Exclusion at the 95% C.L. means that the Higgs hypothesis has at most 5% probability. If no Higgs 
signal has been found, the lower the predicted signal cross section, the more probable is the Higgs 
hypothesis. Hence, for given mass mux and search channel for the Higgs boson, one can determine 
a lower bound Qmm imm ), for cross sections above which the probability of the Higgs hypothesis 
stays below 5%. A model with cross section (T m odel(w#i) > C 7 m i„ (mm ) for the given search chan- 
nel is then said to be excluded at the 95% C.L. Examples of typical (here: LEP) search results are 
shown in Fig. [j] Each of them assumes a particular Higgs decay pattern: (a) SM decay branching 
ratios, (b) 100% decay into bb, and (c) 100% decay into t + t~. As the decay pattern is fixed, the 
limit is given on the Higgs production cross section (Higgsstrahlung) only (and normalised to the 
SM prediction): S 95 (m m ) := a^(m m )/a^(m m ). 

For the SM result in Fig. [j](a), several individual Higgs search topologies have been combined 
in order to increase the sensitivity to a SM Higgs boson. Hence, the result applies only to SM-like 
models, i.e. models where all individual signal cross sections differ from the SM at most by the 
multiplication with one common proportionality factor. Fig. |l|(b) and Fig. [T](c) show results for 
single search topologies, which can be applied to any model which predicts Higgs production via 




Figure 1: Upper bound on the expected (dashed) and observed (solid) cross section ratio S95 for Higgs 
production via Higgsstrahlung, assuming for the Higgs decay: (a) SM branching ratios, (b) 100% decay into 
bb, or (c) 100% decay into t+t~. [§] 
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Higgsstrahlung and has a non-zero BR to bb or t + t , respectively. If for a model which is not 
SM-like a^ z del (m // )BR model (// -> bb)/a^(m m ) > S° 5 bserved (m H ) in Fig. ^b), then the model is 
excluded at the 95% C.L. by this search topology. 

If exclusion results on several individual search topologies are available for a given model, the 
question arises which one to apply in order to determine if it is excluded. Claiming exclusion of 
the model if at least one search topology excludes it certainly does not preserve the 95% C.L. and, 
therefore, should not be called exclusion. The way out is to always consider only one search result 
for excluding a model, but chose the analysis with the highest statistical sensitivity for that. This 
information is encoded in the expected limit (dashed lines in Fig. [[]), the limit which would result 
if only background events were present in the data (obtained via Monte Carlo simulation). 

2. Implementation 

HiggsBounds tests the predictions of models with arbitrary Higgs sectors against exclu- 
sion bounds obtained from LEP/Tevatron Higgs searches. It provides access to all relevant Higgs 
exclusion limits, including expected limits. As input, the program requires the predictions of the 
considered model for: the number of neutral and singly charged Higgs bosons h, , m^, Higgs to- 
tal widths r tot (hj), decay branching ratios BR(/j, —>•...), and production cross section ratios with 
respect to reference values. For a given model scenario, HiggsBounds considers all imple- 
mented LEP and Tevatron analyses, checks whether they can be applied to the model considered, 
and returns whether the scenario is excluded at the 95% C.L. or not. There are three ways to use 
HiggsBounds: a command line version, a subroutine version (written in Fortran 77 and 90) a 
web interface. 

The Tevatron and LEP cross section limits are understood to be applicable to models which 
do not change the signature of the background processes considerably and are usually given for a 
narrow-width Higgs boson. From this and other sources result some limitations on the applicability 
of HiggsBounds which are described in detail in [Q]. 

Now to a crucial bit of the code, namely, how one can make a statement at the 95% C.L. while 
considering several search analyses in parallel. First a definition: we call an "analysis application" 
X, the application of a Higgs search analysis to a particular Higgs boson (or e.g. two Higgs bosons, 
if the analysis involves two of them) of the model under study with particular mass(es). To each 
analysis application X corresponds a signal cross section prediction a(X) for the particular Higgs 
boson on which an upper limit is put by the analysis. For each X, HiggsBounds uses the input 
to calculate the quantity <2modeiP0, which is a(X) up to a normalisation factor. In order to ensure 
the correct statistical interpretation of the results, it is crucial to only consider the experimentally 
observed limit for one particular X. Therefore, HiggsBounds must first determine Xq, which is 
defined as the analysis application with the highest statistical sensitivity for the model point under 
consideration. In order to do this, the program uses the tables of expected experimental limits to 
obtain a quantity Q expec corresponding to each X. The analysis application with the largest value 
of <2modei/<2expec is chosen as Xq. The program then determines the value for <2obs(^o). using 
the appropriate observed limit. If £2modei(^o)2obs(^o) > L HiggsBounds concludes that this 
particular parameter point is excluded at 95 % C.L. 
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In HiggsBounds 2.0.0, 82 Higgs search analyses have been implemented (29 from LEP 
and 53 from the Tevatron). With respect to HiggsBounds 1.2.0 [@], many new types of anal- 
yses have been added and several Tevatron analyses have been replaced by updated ones. The 
following search topologies are implemented in HiggsBounds 2.0.0. (A full list of analyses 
with references is given in the manual 

• LEP neutral Higgs analyses considering the final states: 

hkZ,hk — > bb I TT / anything / invisible I yy I hadrons; bbhk — > bbbb I bbxx (hk CP even or odd) 
; zzhk — > TTTT (hk CP even or odd) ; hkZ,hk — > hjhj,hi — > bb I TT ; /j^/z,-, h^M — > bb 111 ; 
h^hi, hk —> h(hi, hi — > bb I TT ; hkZ,hk — > hjhi,hj — > bb, TT ; h^hj, hk — > bb,h{ — > TT. 

• Tevatron neutral Higgs single topology considering the final states: 

Zhu — > Ubb ; Wht — > Ivbb ; bhk -> 3b jets ; single hk — > WW ; single hk — > TT ; Whk — > 3W ; 
bhk — > bxx ; tthk — > ttbb ; single hk — >■ Zy. 

• Tevatron SM Higgs combined topologies considering the final states (schematic): 

Vh k ^bb + Ef ss (V = W,Z) ; hk+X^WW + X; h k ^WW^ll; h k + X^rr; h k +X^ 
bb + X; Vhk^VVV ^l ± l ± +X (I = e,n,V = W,Z) ; h k +X -^yy+X; h k +X. 

• LEP/Tevatron charged Higgs analyses: 

e + e~ -f H + H 4 jets ; e + e~ -> H + H tvtv ; pp ^tt,t H + b, H + cs (& c.c.) ; 
pp ^tt,t ->■ H + b, H + TV (& c.c). 

3. Applications 

• 4 th Generation Model (FGM): As a toy example, we study a simplified FGM, where we fix 
T(H — > gg) mo dei = 9T(H — > gg)sM for all Higgs masses 1 . Fig. §(A) shows the ratio of the signal 
cross section limit vs. the model prediction for the SM (solid) and the FGM (dashed) obtained with 
HiggsBounds 2.0.0. The SM curve reproduces nicely the LEP and Tevatron exclusions pre- 
sented at the ICHEP' 10. In the FGM, o(pp — > gg — > H — > W + W~) is strongly enhanced compared 
to the SM. Hence, the Tevatron Higgs mass exclusion in this model stretches from 130 GeV to 220 
GeV. Given our simplified approach, this is in very good agreement with [Qj. 

• MSSM Benchmark Scenarios: HiggsBounds is useful to update LEP exclusion plots for 
MSSM benchmark scenarios. Fig. ^(B) shows the good agreement between LEP results for the 
m™ ax scenario [Q] (left panel) and results obtained with HiggsBounds 1.0.0 (right panel) with 
updated MSSM predictions (using FeynHiggs 2.6.4 [J7|]). Verison 1.0.0 already contains all 
relevant LEP search channels for this comparison. Slight differences originate from the fact that 
HiggsBounds relies on LEP data for single topological cross sections, while the full LEP analy- 
sis combines the different search channels. Fig. |3| shows for the MSSM wi™ ax scenario the exclusion 
region (left panel) and the analyses with highest statistical sensitivity in the — tan j8 plane ob- 
tained with HiggsBounds 2.0.0, i.e. including all Tevatron analyses presented at the ICHEP' 10, 
and FeynHiggs 2.7.1 [^] for the model predictions. The gap in the Tevatron exclusion for 
200 GeV <m A < 220 GeV could possibly be filled if the CDF & D0 combined tt analysis [§], 
which has the highest sensitivity for the Tevatron exclusion for HOGeV < < 200GeV, could 
be extended to values m& > 200 GeV. 

'Assuming very heavy 4 th generation quarks, an enhancement factor of 9 is only valid for mu <S 2m t , i.e. towards 
large mu in Fig. 0(A), this overestimates the signal cross section by 20%. 
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Figure 2: (A) Ratio [a x BR] G b s . limit /[<? x BR] mo dei for the most sensitive analysis as a function of the 
Higgs mass niH'. 4th Generation Model versus SM. (B) Excluded region in the to/, — tan/3 plane for the 
MSSM m max scenario. Results from LEP [|| and obtained by using HiggsBounds 1.0.0 are compared. 



• Scalar Sector of the Randall-Sundrum Model: The Randall-Sundrum (RS) Model considers 
spacetime a slice of 5d anti-de-Sitter space with two boundaries, the IR brane (our 4d spacetime) 
and the UV brane [Q], which naturally explains the hierarchy problem. As a consequence of stabil- 
ising the compactification "radius" in the model, an additional scalar, the radion <p , appears in the 
spectrum [10]. Higgs-radion mixing may occur via the interaction = — E, ^—gm&R{gmd) 
where gi n d is the induced 4d metric on the IR brane, R the Ricci scalar, and <I> the Higgs doublet. 
Hence, in general, the physical Higgs Iiq and the radion <po mix to form two mass eigenstates h, (p. 
Like for the Higgs, the radion couplings to massive fermions and gauge bosons are oc mass, but 
e.g. the couplings (pobb and <po 77 are suppressed while (pogg is enhanced with respect to the SM 
Higgs. Although presently many modifications to the original RS Model (RSI) are considered, the 
RS scalar sector is a rather robust prediction, as it is closely linked to the solution of the hierarchy 
problem. Fig. |] shows for RS 1 the excluded region (left panel) and the search channel with highest 
sensitivity (right panel) in the — plane, while the two other free parameters in this model are 
set to A,p = ITeV and = 1/6. This scenario shows slight Higgs-radion mixing, i.e. h behaves 
mainly like the SM Higgs and q> mainly like the unmixed radion. Given the properties of ho and <po, 
this is reflected in several ways by the pattern of LEP/Tevatron exclusions and highest-sensitivity 
channels as a function of m/, and m^. 



4. Summary 

The Higgs search at Tevatron and LEP turn(ed) out many limits on cross sections of individual 
and combined signal topologies. Published individually, it is no easy to make use of all of them. 
HiggsBounds 2.0.0 is a model-independent tool which offers easy access to a wealth of pub- 
lished limits in three ways: command line, subroutines, web interface. It offers a flexible range 
of input formats for the necessary model predictions, including the number of neutral and charged 
Higgs bosons. Applications to a 4 th generation model, the MSSM and the Randall-Sundrum model 
demonstrate its usefulness. The code is publicly available [Q]. 
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■ : pp -> H/A -> r + T~ [CDF '09] 

■ : pp -> K/v4 -> T + T~ [DO note 5740] 

□ : pp — > b h/A — > br + T~ [DO note 6083] 

□ : pp — > h/A — > t + t~ [CDF & do 10] 

■ : pp — > H/A — > T + T~ [CDF & DO '10] 
x: pp — > h/H/A -> T + T~ [CDF & DO '10] 



Figure 3: Excluded region in the tha — tan/3 plane for the MSSM m 
sensitive analysis (right panel) using HiggsBounds 2.0.0. 
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Figure 4: Excluded region in the m/, — m,, plane of the scalar sector of the original Randall-Sundrum model 
(RSI) (left panel) and most sensitive search channel (right panel) using HiggsBounds 2.0.0. 
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